Emergent gravity in two dimensions 
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We explore models with emergent gravity and metric by means of numerical simulations. A par- 
ticular type of two-dimensional non- linear sigma- model is regularized and discretized on a quadratic 
lattice. It is characterized by lattice diffeomorphism invariance which ensures in the continuum limit 
the symmetry of general coordinate transformations. We observe a collective order parameter with 
properties of a metric, showing Minkowski or euclidean signature. The correlation functions of the 
metric reveal an interesting long-distance behavior with power-like decay. This universal critical 
behavior occurs without tuning of parameters and thus constitutes an example of "self-tuned criti- 
cality" for this type of sigma-models. We also find a non- vanishing expectation value of a "zweibein" 
related to the "internal" degrees of freedom of the scalar field, again with long-range correlations. 
The metric is well described as a composite of the zweibein. A scalar condensate breaks euclidean 
rotation symmetry. 



I. INTRODUCTION 



There are many attempts to formulate a quantum field 
^heory for gravity by use of a lattice regularization. Regge- 
Wheeler lattice gravity [lj employs the lenghts of edges of 
+ 7 > simplices as basic degrees of freedom and therefore uses 
Qjiirectly elements of (discrete) geometry. Different geo- 
^ metrical objects are used in other formulations of lattice 
^^gravityJUQ. Lattice spinor gravity Q is an example where 
the basic degrees of freedom are fermions, while no basic 
O^l geometrical objects are introduced. An extension of lattice 
spinor gravity uses in addition to the spinors a geometrical 
field, namely a connection 0. 

t-H We follow here the approach that the metric is obtained 
as the expectation value of a suitable collective field, while 
qq* geometrical quantities are not used as fundamental degrees 
of freedom 0,0] • In this sense gravity and geometry emerge 
{N) from a model of other, non-geometrical fields. In our ap- 
proach, the decisive ingredient is diffeomorphism symmc- 
^Jtry in the continuum limit. In four dimensions, diffeo- 
morphism symmetry entails under rather general circum- 
stances the presence of a massless spin-two particle, and 
Jrj therefore of a metric and the corresponding geometry. We 
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employ a lattice formulation with the property of lattice 
diffeomorphism invariance [8( of the action and functional 
measure. This induces diffeomorphism symmetry whenever 
the model exhibits non-trivial long range physics which al- 
lows to formulate a continuum limit. 

Lattice diffeomorphism invariant models have been pro- 
posed using fermions as fundamental degrees of freedom 
- namely lattice spinor gravity Q. A different approach 
is formulated as a non-linear a- model Q. This latter ap- 
proach offers the important advantage that relatively cheap 
numerical simulations can be used in order to compute ex- 
pectation values and correlation functions of the collective 
metric field. The purpose of the present paper is a first 
numerical study of such lattice diffeomorphism invariant 
non-linear cr-models. 

Our first approach concentrates on two dimensions. 
Sometimes it is said that two-dimensional gravity is triv- 



ial since it does not exhibit a propagating degree of free- 
dom. This statement holds, however, only if the quantum 
effective action is given by an Einstein-Hilbert term pro- 
portional to the curvature scalar R, plus a cosmological 
constant. Indeed, in two dimensions R is a topological 
invariant which cannot provide a kinetic term for the met- 
ric. We emphasize that there is no reason to believe that 
the quantum effective action for a composite metric should 
be purely of the Einstein-Hilbert type. Propagating met- 
ric degrees of freedom become possible for a different form 
of the quantum effective action. We demonstrate this in 
appendix |A"1 with a rather simple possible form of a diffeo- 
morphism symmetric effective action. Similarly, we give in 
appendix [B] an example for a diffeomorphism symmetric 
effective action for a propagating zweibein. The effective 
action for the model investigated in the present paper turns 
out to be more complicated than these simple examples. 

In this paper we investigate a non-linear er-model with 
two complex scalar fields cpi, with "flavor" index i — 1,2. 
The constraint 
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(1) 



is compatible with an 50(4)-flavor symmetry acting on the 
four real fields ipi,R,ipij,(pi — <Pi,R+i<Pi,i- We consider the 
real action (x M = (a; ,^ 1 )) 



S = (3 



with e 01 



d 2 xe^ 



(2) 

= I, e 00 = e 11 = 0. For continuous space 
this action is invariant under general coordinate transfor- 
mations. No metric is introduced a priori - diffeomorphism 
symmetry is realized by the particular contraction of two 
derivatives with the e-tensor. 

We discretise the action )| on a two dimensional 
quadratic lattice by using lattice derivatives and cell aver- 
ages as explained in detail in section ITT] We perform a nu- 
merical study of this model by Monte Carlo technics. The 
continuum diffeomorphism symmetry of the action implies 
diffeomorphism symmetry of the discretised lattice action. 
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However, we stress that the effective action describing the 
system in the continuum limit is not of the simple form 
Our main findings are the following: 

(i) We identify collective fields that have the transforma- 
tion properties of a metric . They acquire indeed non- 
vanishing expectation values. There are several can- 
didates for collective metric fields. For the vacuum we 
find non- vanishing flat metrics both with Minkowski 
and euclidean signature. The metric turns out to de- 
scribe flat space independently of the parameters of 
the model, which raises interesting questions concern- 
ing a self- adjustment of the effective two-dimensional 
cosmological constant to zero. 

(ii) The correlations of the metric fluctuations are long 
range. They typically show a powerlike decay ~ r~ Q , 
with a close to two. 

(iii) The geometry differs from flat space space if sources 
coresponding to an energy momentum tensor are in- 
troduced. In the linear approxiamtion the pertur- 
bation of the metric in response to a source is de- 
termined by the correlation function. A static point 
source, which mimicks a static massive object, leads 
to a "Newtonian" potential that decays with the in- 
verse of the spatial distance. This is similar to four 
dimensional gravity, but quite different from naive 
dimension estimates. 

(iv) The critical behavior associated to the powerlaw for 
the correlation functions occurs independently of the 
detailed values of the parameter /3 characterizing the 
model. The non-linear cr-model © is an example for 
self-tuned criticality. 

(v) We also identify a collective zweibein e™, where the 
"Lorentz index" m is associated to the flavor degrees 
of the scalar field. To a good approximation, some of 
the metric candidates can be described as the usual 
bilinear of the zweibein. 

(vi) The correlation functions of the zweibein decay with 
a power of the distance, similar to the metric. 

(vii) We identify the exact ground state of the model for 
(3 — > oo. It consists of "stripe-configurations" for the 
scalar field. 

(viii) The stripes persist for j3 > f3 c ,f3 c rj 4.6. They 
also characterize the vacuum state of a correspond- 
ing continuum theory. The scalar order parameters 
describing the stripes are responsible for a sponta- 
neous breaking of euclidean rotation symmetry. For 
this reason we encounter, in case of euclidean sig- 
nature, an unusual version of d = 2 gravity. While 
diffeomorphism symmetry is expected to be realized 
in the continuum limit this does not hold for rotation 
symmetry, due to the presence of preferred axes. 

(ix) At /3 C we find a first order transition. The disordered 
phase for j3 < j3 c shows no stripes and no expectation 



value for the zweibein. The correlations are short 
range in this phase. 

Our paper is organized as follows. In section [TT] we de- 
scribe the discretisation of our model on a lattice. We de- 
fine scalar order parameters and collective metric fields in 
sect. IIIII In section HVl we show that the metric correlation 
functions have a long range, power law behaviour. The 
stripe configurations that dominate the system for large 
values of j3 are discussed in sect. [Vj We find the exact 
"ground state" of the system for j3 — > oo. It is a stripe 
configuration for the scalar fields. In section IVTl we define 
and describe the zweibein as a bilinear in the scalar fields. 
The "internal index" or "Lorentz-index" of the zweibein 
is related to the flavor structure for the scalars. In sec- 
tion lVIll we discuss the presence of an approximate Lorentz 
symmetry of our non-linear cr-model. While the action ^ 
is Lorentz symmetric, the constraint (JlJ violates Lorentz 
symmetry. We describe the metric and zweibein correla- 
tion functions in sections fVIIII and llXl respectively. Finally, 
we conclude in section [X] 

The results presented in this paper are all obtained by 
numerical simulations. In parallel, we present in a series 
of appendices our analytical investigations of the possible 
form of effective action for two-dimensional gravity. These 
analytical considerations all apply to the continuum limit 
and exploit the symmetries of our model. We proceed on 
various levels: App. IA1 discusses the metric as a unique de- 
gree of freedom, while App. IBluses the zweibein. In App. [D] 
we investigate the effective action for the scalar fields, and 
App. lEl combines scalars and the zweibein. Only the effec- 
tive actions in App. [TJ1 and lEl describe parts of our numer- 
ical findings in a realistic way. Finally, App. [C] discusses 
a possible connection between the correlation functions for 
the zweibein and the one for the scalars. 



II. LATTICE ACTION 

We regularize our model on a lattice, x° = z°A, x 1 = 
z 1 A, with z M integers such that the sum z° + z 1 is odd. 
This corresponds to a square lattice with lattice distance 
\/2A and nearest neighbors in the diagonal directions ~ 
x° ± x 1 . For the corresponding "diagonal lattice vectors" 
the z-coordinates are 

So = (1,1) , Si = (-1,1). (3) 

We define cells located on the sites of the dual lattice at 
x fj- — y^A, with y M integer and y° + y 1 even. Each cell 
consists of four lattice points with lattice coordinates 5^ = 
y^ ± (vv)^, with two unit vectors v v obeying (v u )^ = 6%. 
The lattice derivative at y is given by 

d^W) = ^ [<p(ii + - <p(v - > ( 4 ) 

and cell averages obey 

= 2[v(v+ v o)+v(y-Vo)+<p(y+Vi) + (p(y-Vij\. (5) 
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The lattice action replaces in eq. ([2]) the derivatives by 
lattice derivatives, the fields without derivatives by the cell 
averages, and J d 2 x —> J2y V(y) with cell volume V(y) = 
2A 2 . The distance between neighboring cells is \/2A. 

Most important for our context, the action is lattice dif- 
feomorphism invariant. Indeed we can change the posi- 
tioning of the lattice points x^(z) from the regular lat- 
tice x^(z) = Az p to an arbitrary neighboring lattice 
x^(z) = Az^ 1 + ^(z). Here the cartesian coordinates x^ 
parametrize a manifold which is some region in R 2 . For 
general positions of the lattice points on the manifold the 
lattice derivatives read [8j 

do<fi(y) = ^^c^(x 1 {y + vi)-x 1 {y-vi)) 

(<p(y + vo) - <p(y - v )) 

- {^(y + v ) -x 1 ^)- v )) (tp(y + vi) - tp(y - ui))}, 

9i<p(y) = 7j^^{( x °(y + v o)-x°(y-vo)) 

(ip(y + v x ) -tp{y-vx)) (6) 

- (x°(y + vi) - x°(y - vi)((p(y + v ) - ip(y - «o))j, 

and the cell volume becomes 

V(y) = ^e IJlV (x' 1 {y+Vo)-x' i {y-VQ))(x v '{y+v x )-x v '{y-v x )). 

(7) 

The cell averages ([5]) remain the same. The change of 
each term in the action due to the change of the lattice 
derivatives is precisely canceled by the change of the vol- 
ume factor in J d 2 x = ^yV{y). The expression of the 
lattice action in terms of lattice derivatives and cell aver- 
ages is therefore independent of the positioning of the lat- 
tice points. This crucial property is due to the particular 
contraction of the lattice derivatives with the e-tensor. It 
guarantees the standard diffeomorphism symmetry of the 
quantum effective action in the continuum limit |8j. 

We will impose periodic boundary conditions on a square 
lattice defined by the Eo and E\ vectors. We have checked 
by performing calculations on a square lattice defined by 
the Xo and x\ axes, that the boundary conditions do not 
play an important role in determining the phase structure 
of the theory. The functional integral 

Z = J Vipe~ s (8) 

involves for every lattice point z the standard 50(4) in- 
variant measure on the sphere S 3 in the field space. The 
functional integral is finite for a finite number of lattice 
points. Our model is therefore mathematically well defined 
and a candidate for regularized quantum gravity. 

We work with a fixed positioning of the lattice points 
on the regular quadratic lattice as described above. For 
general relativity, this corresponds to a fixed choice of co- 
ordinates. 

We obtain the numerical results presented in this paper 
by Monte Carlo simulations using the standard Metropolis 



algorithm on 64 2 or 256 2 lattices. We usually start the 
thermalization process in the stripe phase (see sect. N} . 
We have checked that the system eventually ends up in 
the stripe phase even if we start the thermalization in the 
disordered phase (for the lattice sizes employed here). In 
this case the thermalization typically takes much longer, 
as the system first breaks up into domains with different 
orientations of stripes, which slowly equilibrate into one 
coherent domain. 



III. ORDER PARAMETERS AND METRICS 

The continuum action ([2]) is invariant under translations 
and rotations, as well as with respect to a parity-type dis- 
crete transformation P: x 1 — > —x , ipi f-> ip2 or time rever- 
sal T : x° — > —x°, ipi f-> <^2- Similarly, S does not change 
under diagonal reflections D± : x° f-> ±ar, (fi -H> (f2- 
Another discrete symmetry is charge conjugation which 
is realized by complex conjugation, C : ifi — > ip*. Fur- 
thermore, the action conserves a global continuous flavor 
symmetry with abelian gauge group U(l) x C/(l), corre- 
sponding to separate phase rotations for ipi and f>2- Since 
the action changes sign under tpi f-> if2 we can restrict the 
discussion to positive (3. We also observe that S changes 
sign if only one of the fields is replaced by its complex 
conjugate, say ipi — > (p\",<p2 — > ^2- For the reflections 
P, T, D± we can therefore replace the accompanying reflec- 
tion ifi -H- if2 by ip\ <H- Lp\. Combined flavor reflections, as 
fi f2, <y52 - > ¥>i, leave the action invariant. 

The lattice action is invariant under diagonal transla- 
tions of V2A, and therefore also under translations in 
the x° or x 1 directions by 2 A. It is preserved by ir/2- 
rotations and shares the same discrete symmetries as dis- 
cussed before for continuous space, as well as the continu- 
ous U(l) x U(l) flavor symmetry. 

We find that for large enough j3 several of the discrete 
symmetries are spontaneously broken. Indeed, the char- 
acteristic configurations of the scalar fields change qual- 
itatively at (3 C w 4.6, as demonstrated in Fig. [1] For 
fi > f3 c one observes an order in stripes which disappears 
for /3 < f3 c . For the "stripe phase" (f3 > f3 c ) we can define 
order parameters sf by defining "supercells" with sixteen 
lattice points. For sf the field values ipi within a given su- 
percell are summed with phases according to the left part 
of Fig. [21 while for s~ we sum with the complex conjugate 
of the phase factors, 

z 

s r = i2E eV(5 Vi(2)- (9) 

z 

(More precisely, the phases exp(ir/(z)) depend on the coor- 
dinates within the supercell, as defined by the left panel of 
Fig. [2] and L is the number of lattice points along one (di- 
agonal) direction.) We will show in sect. |V]that there are 
four different classes of stripe configurations which can not 
be rotated into each other using the internal C/(l) x U(l) 
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symmetry. Each of the four order parameters sf with 
i = 1,2 and a = ± signals the realization of one of the 
four equivalent classes of stripe configurations. 

We show a characteristic order parameter as a function 
of (3 in Fig. [3] In the presence of such order the lattice 
translation symmetry is partially broken - only transla- 
tions by 4a/2A in the diagonal directions or by 8A in the 
x° or x 1 directions leave the equilibrium state invariant. 
Also the discrete reflection symmetries and the symmetry 
of 7r/2-lattice rotations are broken spontaneously by the di- 
rection of the stripes. The symmetry of the continuous fla- 
vor transformations gets broken spontaneously as well and 
one may expect Goldstone-type excitations. All these sym- 
metries are preserved in the disordered phase for (3 < f3 c . 
The discontinuity in the scalar order parameter visible in 
Fig. [3] indicates a first order phase transition. This will be 
confirmed by jump in other expectation values. 

For (3 > f3 c we will see in sect. [V] that lattice transla- 
tions, reflections and rotations leave the equilibrium state 
invariant if they are combined with appropriate phase (or 
flavor) rotations. It is therefore perhaps more appropriate 
to associate the stripe phase with a spontaneous breaking 
of the U(l) x {/(l)-flavor symmetry. (We discuss compati- 
bility with the Mermin- Wagner theorem at the end of sect. 
IVD Observables which are invariant under flavor rotations 
will be invariant under translations by Eq or E±. 

We are interested in geometry and therefore look for ex- 
pectation values of observables that can play the role of 
a metric. In the continuum limit such observables should 
transform as second rank symmetric tensors. We will find 
several natural candidates for composite metric observ- 
ables, with euclidean or Minkowski signature of the ex- 
pectation values. They are typically invariant under flavor 
rotations. 

Metric tensors can be constructed from derivatives of cpi. 
Many possibilities exist for constructing objects that trans- 
form as a symmetric second rank tensors. In the contin- 
uum limit ip transforms under diffeomorphisms as a scalar 
and d^ip as a covariant vector. The symmetric product 
of two vectors transforms therefore as symmetric second 
rank covariant tensor, which is precisely the transformation 
property of the metric. The lattice analogon replaces the 
derivatives by lattice derivatives ((4]) and the fields without 
derivatives by the cell averages ([5]). We discuss here four 
examples that are all invariant under the flavor symmetry 
of continuous phase rotations of ipi and if 2- 

Our first candidate for a metric reads 

g$(y) = 2A 2 Re(Y,d^*(y)d^(y)), (10) 

i 

with lattice derivatives d^<f(y) given by eq. (|4]). On a L x L 
lattice we find that the expectation value is invariant under 
lattice translations, 

9$(y) = (g$m = NW(f3)5,„, (ii) 

with N^ 2 ' (f3) shown in Fig. 2J The corresponding geometry 
is euclidean flat space for all (3. (The normalization factor 
((3) can be absorbed by a rescaling of the metric or the 



coordinates.) We observe a discontinuity of N( 2 \(3) at (3 C . 
The particular ground state g^ v ~ 5 M „ preserves the lattice 
rotations and reflections. 
A second candidate is 

gtfP = ^ 2 Re{iipld^ x )Re{iipld v ^) + {<pi <-► fflO- (12) 

For this "Minkowski metric" we find that the expectation 
value 

C ) = (~9$ ) )=NM(f3)r ) » u (13) 

has Minkowski signature r]^ — diag(— 1, 1). Again, 
((3) is plotted in Fig. 2] and shows a discontinuity 
at f3 c . The expectation value (fT3"]) is the same for all lattice 
sites. It is invariant under (euclidean) rotations by 7r/2 and 
under the discrete symmetries P and T. 
A third metric can be written as 

g$ D > = %A 2 {Re{if\d li f X )Re(if\d v f 1 ) - {<fx <p 2 )}. 

(14) 

The expectation value gffi = (9^°^) is found as 

9g m =g£ D) = o,&' D >=Ni»*w. (15) 

Again, this corresponds to a Minkowski signature since 
det{gif D) ) = -{N^ MD \I3)) 2 . In the continuum, the met- 
ric (jT5"j) obtains from the metric (fT3")l by a euclidean rota- 
tion of the coordinate axes by 7r/4. For the metric (|15[) 
the time- and space- axes are given by the diagonal axes in 
the directions of Eq = (1,1) and E\ = (—1,1). Finally, we 
investigate a fourth metric 

= ^Reinpld^Re^ld^) + {ip x O <p 2 ). (16) 

The expectation value corresponds again to euclidean flat 
space 

It is remarkable that, according to the precise definition 
of the metric observable, one finds expectation values that 
correspond to flat space, either with a euclidean signature 
or a Minkowski signature. The signature apparently de- 
pends on the flavor structure of the metric observable. We 
will gain later a better understanding of this issue. The 
dependence of the metric observable on j3 shown in Fig. 
13] conforms the discontinuity at the critical value /? c that 
signals the presence of a first order transition. 



IV. LONG RANGE METRIC CORRELATIONS 

In the stripe phase for f3 > f3 c we observe the presence of 
long range correlations for fluctuations of the metric field. 
They decay with a power law. We define the fluctuations 

h®=g$-g$=g$-<3$) (18) 
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FIG. 1: The real part of <pi for a typical configuration at /3 = 
4 < fi c (top) and j3 = 5 > /3 C (bottom) on a 64 2 lattice. Units 
along the axes are y/2A. 





FIG. 2: Phases of fields in a supercell, 1 for +, —1 for — , i for 
• — i for O- The small dot in the center is the location of the 
supercell. 



for the different metric fields dTUj) , (12]), (114]), (JT6J) , labeled 
by (^4). As an example, we consider the correlation func- 
tion 

Gf{r) = \ ((h®(x) - h$( x) )(hg\y) - hg)(y)) 



(19) 

for r — \x — y\. In Fig. [5] we plot G d (r) on the diagonal 
axis ~ Eq for three values of j3. For (3 > /3 C we clearly see a 
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FIG. 3: The order parameter \sf\' 2 as a function of /3. The 
thermalization of the system was started from the striped phase 
signalled by . 
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FIG. 4: Proportionality factor N(/3) of the metric expectation 
values defined in eqs. (|10p . (|12[) . (I14p and (|16|l . as a function of 



power law decay. A fit for the range 5A < r < 50A yields 
Gf{r) =c d r~ a * (20) 



with 



a d = 2.07 ±0.05 
a d = 1.99 ±0.03 



c d = 0.031 for /3 = 4.7 

c d = 0.0069 for /3 = 20. (21) 



It is remarkable that the powerlaw decay of the corre- 
lation function occurs for all values (3 > f3 c . No tuning of 
a parameter, as common for critical behavior at a second 
phase transition, is necessary. We have excluded from the 
fit the points with r < 5A since lattice details are relevant 
in this range - a possible universal behavior will only be 
found at distances r sufficiently large compared to the lat- 
tice distance. For r > 50A the correlation becomes very 
small and the uncertainties large. Also finite volume effects 
may start to play a role. 

On the other hand, we observe in the symmetric phase 
for (3 < f3 c a much faster decay of the metric correlation. 
There are only few points before the correlation gets very 
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FIG. 5: Euclidean metric correlator G^ 2 (r), defined in eq. (|I9[) . 
as a function of distance along the Eq axis, for different cou- 
plings /3. 



observed long range correlations for the metric observables. 
Further aspects of the correlation function for the metric 
will be discussed in sect. IVIIIl 



V. STRIPES 

For large values of f3 the leading configurations are 
stripes. A characteristic stripe configuration reads 

MS) = ^jf aAl) > = 7f M2(5) ' ( 23 ) 



with phases 



oti (y + (2m 1 - l)v + riiEi) ■■ 
a 2 (y + (2m 2 - l)vi + n 2 E ) = - 



TO27T 



(24) 



small and the uncertainties large. The decay is compatible 
with an exponential decay or a power law with large nega- 
tive exponent. Universal long distance behavior seems not 
to be realized in the symmetric phase, even for /3 close to 
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Here m$ and rij are arbitrary integers and Eq, E\ are given 
by eq. ((HI)- For ip.\ the stripe is in the diagonal ^-direction, 
with tpi(z + nEi) = </?i(5), while Lpi it is orthogonal to it 
in the diagonal i?o-direction, with ipi(z + uEq) = Lpi(z). 
A typical realistic stripe configuration with fluctuations is 
shown in the right part of Fig. Q] While tpi is invariant 
under i?i-translations, it is invariant under translations in 
the i^o-direction only by four units, <pi(z + AuEq) — ipi(z). 
Similarly, ip% is invariant under -E^-translations by four 
units, ipziz + AnEi) = (p 2 (z). As a whole, the stripe config- 
uration (|23|) . (|24l) is therefore invariant under translations 
in the E - and i^i-directions by four units. This results in 
translation invariance in the a; - and a^-directions by 8A. 

A given stripe configuration is not invariant under the 
continuous phase rotations of the U(l) X [/(l)-flavor sym- 
metry. The phase changes can be used, however, to estab- 
lish a symmetry of translations in the Eq- and ^-directions 
by one unit combined with an appropriate phase rotation. 
Indeed, the combined diagonal translations, 

t : <p'x(z) = -i<fi{z — Eq) , tp' 2 {z) = (p 2 {z — Eq), (25) 



FIG. 6: Minkowski metric correlator G D (r), defined in eq. 
(|22p as a function of distance along the Eq axis, for different 
couplings p. 



The correlation function for the fluctuations ' of the 
Minkowski metric also shows a long range behavior for j3 > 
j3 c . In Fig. [S] we display G^\r), defined by 

<#°(r) = \ ((h f\ X ) + h[Y\ X )) (hif\y) + h[f\y))) 

(22) 

for r — \x — y\ again on the diagonal axis ~ Eq. This 
correlation function shows a power law decay similar to the 

(2) 

euclidean metric gj^J , with a decay exponent close to a = 2. 
In the following we will concentrate on the stripe phase for 
/3 > f3 c . Our aim is an understanding of the nature of the 



and 

h : <p'i(z) = <pi(z — Ei) , ip 2 (z) — —iip 2 (z — Ei), (26) 

leave the stripe configuration (|2"5|) . ([M]) invariant. By 
virtue of these combined translations it is often sufficient to 
evaluate the value that a cell-observable takes for the stripe 
configuration only for the cell at y = 0. The value for the 
neighboring cell configuration at y — Eq obtains from the 
value of the observable for the cell y — by multiplying 
each factor ipi by a phase factor —i. Similarly, the value 
in the neighboring cell at y = E\ is found by multiplying 
each factor tp 2 by a phase —i. This can be continued for 
all cells. In particular, all cell-observables that are invari- 
ant under U(l) x L/(l) flavor transformations take for the 
stripe configuration the same value in each cell. 

We first evaluate for the stripe configuration the value 
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of the cell action 



S = 2A 2 /3 (27) 



This cell observable is independent of flavor-phase rota- 
tions and takes therefore the same values in all cells. For 
y = one finds for the stripe configuration ((23]), (1241 the 
following values 



1 . . 1 . 

(1 - I) , tp 2 = — 7=(1 ~ l), 



2V2 



2V2 



dotpi = - n / - A (1 + i) , = - 



2v^A ' 2V2A 

' (1 + i) , &^ = -_L-(i+i). (28) 



2\/2A 



2\/2A 



This yields 



<pld <fii = - 



1 

4A 

i 

4A 



(plditpi = - 
<P2<h<P2 = - 



1 

4A' 

i 

4A' 



and therefore 



2 A 2 ' 



(29) 



(30) 



We will see below that the stripe configuration (|2"3")l . 
minimizes the action ([2]). (The action per lattice point 
takes the value — j3.) The ground state for j3 — > 00 is there- 
fore given by the stripe (|23|) . (|24l) or one of the equivalent 
configurations that we will discuss below. 

The metric bilinear (flOl) is again a cell observable that 
does not depend on the flavor phases and therefore takes 
the same value for all cells. For the stripe configuration 
one finds from eq. 



(31) 



This agrees with the asymptotic value N( 2 \/3 — > 00) = 1 
that can be seen in Fig. SJ 

We next investigate the discrete lattice reflections which 
leave the stripe configuration invariant. A time reflection 
at the axis z° = M leaves the stripe invariant if it is ac- 
companied by ipi 4-> if2 and an appropriate global phase 
change 



T, 



M 



where 



- M, S 1 ) = e iaM <p 2 {- (z - M), h), 
ip' 2 {z°-M, h) =e~ iaM <pi(- (zo- M),zi), 



Um 



Mtt 



(32) 



(33) 



Similarly, a parity type invariance of the stripe configura- 
tion obtains by combining the reflection at the axis z 1 = N 
with an appropriate phase shift 



N 



N) = e la »<p* 2 (z°,~(z 1 -N)), 
z'-N) = e la »<pl(z°,-(z 1 -N)), 
(N+l)n 



(34) 



Since the symmetry Ppj in eq. (|34[) also involves the com- 
plex conjugation of the fields ipi it may be associated with 
a CP transformation. Both symmetries P/v and Tjv are 
symmetries of the action. For a ground state preserving 
these symmetries an observable and its associated reflected 
observable must have the same expectation value. 

In particular, we may consider observables that are in- 
variant under the transformations ipi o if2,<Pi ^> as 
well as under flavor phase rotations. Then the expectation 
values of observables for which only the coordinates are 
reflected must be the same as the ones for the original ob- 
servables. The metric observable (JTUJ) is of this type. The 

(2) (2) 
symmetries Tm and P/v imply g^' (y) — since is odd 

under the reflection of one coordinate. 

For the diagonal reflections D± we restrict the discussion 
to reflections at axes through the origin tje — (0,0). Re- 
flections at shifted axes can be obtained by a combination 
of those "basic reflections" with the translations to and t\ 
given by eqs. (|2"5j) . ([2"B)) . (Also the reflections Tm and P/v 
can be related in this way to suitable basis reflections at 
axes through the origin.) The stripe (|2"3"j) . (|2"4l is indeed 
invariant under the diagonal reflections 

D+ ■ <fi(y) = P\{D+y) , <p' 2 (y) = -iip* 2 (D+y), 

D- ■ v'M = -i<pl(D-y) , v' 2 (y) = MD-y), (35) 



with 



(36) 



The symmetries D + and D_ are symmetries of the ac- 
tion. If they are preserved by the ground state they imply 



g ( n>(y) 



(2) 

Now g y Q1 ' 



reflections, while the difference d = 4(<7oo 



is even under the diagonal 



■gfj) is odd. 



Rotations by it/ 2 can be obtained by combining reflec- 
tions. For example, the combination T D + amounts to a 
rotation around the origin with an angle —tt/2, 



TD + y=(y\-f) 



(37) 



With respect to the ir/2- rotations both d and goi are odd, 
while h — goo + gn is even. 

Symmetries of the action that do not leave the stripe 
configuration (|23|). (|24j) invariant lead to equivalent stripes. 
Equivalent stripes differ from the stripe (|2U)) . (|2"41 in po- 
sition, orientation and phases, while they have the same 
value of the action. (If the ground state is given by one 
particular stripe configuration the symmetry transforma- 
tions leading to equivalent, but not identical, stripes are 
spontaneously broken.) As an example, the action of the 



8 



stripe is not modified if tpi and (p 2 are multiplied by global 
phases e 1 ^ 1 and e 1 ^ 2 , respectively. Similarly, the stripes 
can be displaced by an arbitrary number of units Eq or 
Ei, or they can be rotated by tt/2. Equivalent stripes also 
obtain from reflections P, T, D± if those are accompanied 
by ipi <-> ip2 or ip i f> pj. Pure coordinate reflections of 
stripes do not lead to equivalent stripes, however. For a 
pure coordinate reflection of the stripe (|23p , (j2~4")l the action 
becomes positive C(y) — 1/(2A 2 ) - such reflected stripes 
correspond to a maximum rather than a minimum of the 
action. Among the symmetries of the action that lead to 
equivalent (but not identical ) stripes is the flavor rotation 



Rf : (fx -> (p 2 , (f 2 
and the charge conjugation: 



R t 



(38) 



(39) 



They define equivalence classes of stripes that cannot be 
rotated into each other by phase rotations. While R c maps 
<-> in eq. ([9]), the maps between sf and sf can be 
achieved by Rf. 

For P — > oo the ground state is indeed given by a par- 
ticular stripe and the symmetries of flavor rotations are 
spontaneous broken. This issue is more delicate for finite 
values of (3. Now the Mermin- Wagner theorem [T^] forbids 
any spontaneous breaking of a continuous symmetry in two 
dimensions. Nevertheless, for all practical purposes the 
system behaves for large enough /? as if the flavor sym- 
metries were spontaneously broken. This issue is similar 
to the Kosterlitz-Thouless[llj transition where in the low 
temperature phase a mode with all the properties of a Gold- 
stone boson exists. For the Kosterlitz-Thouless transition 
the general aspects are well understood by applying the 
functional renormalization grou p to linear and non-linear 
cr-models in two dimensions |12j | . The Mermin- Wagner the- 
orem looses its practical applicability. While it remains 
formally valid in the infinite volume limit, the effects of 
spontaneous symmetry breaking occur for a macroscopic 
system with arbitrary (but finite) size. We expect an anal- 
ogous situation for our model, with additional complexity 
due to the lack of rotation or Lorentz-symmetry in the con- 
tinuum limit. 

The limit (3 — > oo projects on the "ground state" for 
which the action takes its minimum value. For our model 
the ground state can be solved exactly - it is the stripe 
configuration (1231) . (j2"4"]l or an equivalent stripe. In order 
to show this we first note that any configuration for which 
C(y) in cq. (|27l) takes its minimum value for every cell 
must be a minimum of the action. We next show that this 
is the case for the stripe. 

Using in eq. (|28j) the definitions (@J|, ©, the cell ac- 
tion 2A 2 £(y) can be written as a sum of terms <p* a (pi3 of 
the eight complex fields ip a in the cell. We want to show 
2A 2 £(y) > — 1. On everyone of the four points in the 
cell the condition ([lj must be obeyed, leaving us with 12 
angles. The flavor symmetry ensures that C(y) does not 
depend on the overall phases of (pi and ip2, such that 10 
angles remain. Finding the minimum of a function of ten 



angles by analytical means is rather involved. We therefore 
have proceeded to a numerical evaluation of 2A 2 £(y) for 
random choices of the angles. For 10 10 configurations we 
have found no value smaller than —1 within the numerical 
accuracy. Together with the value ([30]) for the stripe con- 
figuration we consider this as sufficient evidence that the 
stripes (|2"3"|) . (|24l) minimize the action. 



VI. ZWEIBEIN 



It is possible to express the action ([2]) in terms of 
"zweibein" e M m . We first introduce 



= y/2iA(ftdr<p2 - fadnipZ), (40) 



such that 



2AJ , Sx^E^K 1 



4A 



2A 2 



d 2 xE, 



with 



E = det(£ p 



(41) 



(42) 



In the continuum limit E^ m transforms as a covariant vec- 



tor. The form (|41|) makes the diffcomorphism symmetry of 
the continuum theory particularly apparent. 

It is instructive to evaluate E^ m for the stripe configu- 
ration p5]). (f24|) . One obtains from eq. ([29]) 



p m _ 1 ( 1 , -1 



(43) 



(With E = 1 this reproduces C(y) = -1/(2A 2 ) in eq. 
([3D]).) We may bring the zweibein to a diagonal form by a 
suitable rotation 



For E^ m of the form (|43|) one needs 



R — — t= ( i ' 1 ) i R^ R — 1, 
v/2 V- 1 . 1 



(44) 



(45) 



such that the stripe configuration yields e M m = <5™. We 
will use this convention and define 



7, 



1 



V2 



(E^° - E^ 1 ) = 2ARe(ii P * 1 d^ 1 ) - (^ o tp 2 ) 



eV = ^(^° + E^ 1 ) = 2ARe(ipld^i) + (<pi O <^ 2 ). 



(46) 
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The action retains the form of a determinant of the 
zweibcin 



We have evaluated the expectation value of the zweibein 
and find in the striped phase for /3 > /3 C , 



S 



2A 2 



d 2 xe , e = det(e™)- 



(47) 



A nonvanishing expectation value of the zweibein indi- 
cates the spontaneous breaking of symmetries of the action. 
This is closely related to spontaneous symmetry breaking 
by a given stripe configuration. We can use such symme- 
tries in order to transform among equivalent zweibeins. For 
example, the flavor rotation (|38p induces 



R 



e,,, 1 -> 



-e 



(48) 



(This amounts to a 7r/2-rotation among the flavor indices 
of the zweibein. Comparing with eq. (|45[) we note Rf = 
—R 2 .) Similarly, charge conjugation (f39|) induces 



Rr 



(49) 



(53) 



This holds if the striped phase is charaterized by the con- 
figuration (|23|) . (j24l) . The expectation values in the equiv- 
alent striped phases are given by transformations similar 
to eqs. (|48p and (|4^|) . As usual, a metric can be defined as 
a quadratic expression in the zweibein 



JMe) 



*W = ((tf< e >G8))V- (54) 



It differs from the Minkowski metric (fj^if ' ) , cf. eq. (1521 , by 
the contribution of local zweibein fluctuations. In Fig. [7] we 
display (N^(f3)) 2 together with A< M )(/3). The difference 

is small and we can approximate by gffi with good 
accuracy. 



such that zweibeins with opposite sign are equivalent. 
As usual, the zweibein can be used to define a metric 



(50) 



with fjmn = fjnm- We define the "euclidean metric" g^v 
by choosing f] mn = S mn , 



e °i ° + e x e 1 



(51) 



This coincides with eq. (|16]) . For the stripe configura- 
tion (|23|) . (p4|) one finds = 8^ u . Similar to gffi in 
eq. (fT0|) this metric observable is independent of the flavor 
phases, and remains unchanged for f> ^ or ^ f> 
In contrast to eq. (1101) it involves four powers of fields ipi, 
however. 

Another possibility uses f/ mn = rj mn = diag(—l, 1). The 
metric corresponding to eq. (jf 21) , 
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FIG. 7: Expectation values of the zweibein and metrics, as 
expressed by the proprtionality factors (N ( - e ' > ) 2 ,iV (M) andiV (E) , 
as a function of f3. 



5 (M) _ _~ 0g 



e 1 e 1 



(52) 



takes for the stripe configuration the values g$P — Tj^ u . 
For large enough j3 one therefore expects that the expec- 
tation value (gjfv ) has a Minkowski signature, {g\^P) = 

N^-^ifyVnv We observe that gffi' is again invariant with 
respect to phase changes of ipi or cp2 and to the inter- 
change (fx f-> <^2- It changes sign, however, under the 
exchange ip% O </?*. In consequence, the diagonal reflec- 
tions -D± change g^ 1 ^ ■<-» — ffii while g^i * s °dd- Thus 

„^ (M) . 



Soo° + Sii is even and g$p + g[\ I] odd. 
This additional minus sign extends to 7r/2-rotations under 

which -g$P + g[\ I] and g^ 4) are even and gffl + g[\ I] 
is odd. The expectation values of the metric (|I2p or (|5^|) 
singles out a time direction. This direction depends on 
the stripe orientation. The transformation (|48|) switches 
between time and space directions. Thus the difference be- 
tween space and time occurs as an effect of spontaneous 
symmetry breaking. 



We can go one step further, and try to see if the zweibein 
expectation value is well described by the disconnected part 
of the scalar correlator, that is 



,(¥>) 1 - 



2ARe(i(cp* 1 )(d fl <pi)) 
2ARe(i{<pl){d M tpi)) 



(epi <-> (p 2 ), 



(55) 



Here (ifi) = (<p(y)) denote the averages of the cell aver- 
aged scalar fields, which are space dependent in the broken 
phase, and (d^^y)) are averages of lattice derivatives. 
This quantity has the same structure as e™, 



(56) 



In Fig. we compare AT( e )(/3) to A( e ^(/3), and note that 
the disconnected zweibein describes the full zweibein to a 
reasonable approximation, and the quality of the agree- 
ment gets better at larger (3. 
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FIG. 8: Proportionality factors iV (e) (/3) and N (ev \f3), defined 
in eqs. (|53p and (|56[l . as a function of f3. 



In the continuum limit the presence of order parameters 
e^™ singles out particular directions since both e M ° and 
e^ 1 transform as a vector. A rotation of the axes away 
from the x°- and cc 1 -direction does not leave the zweibein 
diagonal, such that the diagonal form (|53)l singles out a 
particular coordinate frame. The same holds for Lorentz- 
transformations. In the usual vielbein formulation of ge- 
ometry Q the lack of rotation - or Lorentz- symmetry for 
any given value of the vielbein is compensated by an ac- 
companying transformation of the Lorentz frame. A simul- 
taneous rotation between the vectors e M ° and e^ 1 , together 
with a coordinate rotation, leaves the "ground state viel- 
bein" of the type (|53|) invariant. This guarantees rotation 
symmetry for flat euclidean space. A similar property with 
respect to Lorentz symmetry ensures Lorentz symmetry for 
Minkowski space. 

In the Cartan-formulation the euclidean metric g^ v = 
e l / n e v n 8 mn is obtained from the vielbein by contraction 
with the invariant tensor S rnn , such that the metric is rota- 
tion invariant. Then the metric correlation functions trans- 
form under rotations covariantly as dictated by their tensor 
properties, provided that there is no other source of rota- 
tion symmetry breaking. Similar statements hold for the 

Lorentz covariance of the metric (|54p and its corre- 

lation functions since r\ mn is a Lorentz invariant. 

The invariance of the zweibein under combined space 
and internal rotations is, however, not sufficient in order 
to guarantee a rotation invariant setting. In our model 
the expectation value of scalar fields in the form of stripes 
breaks the rotation symmetry without the possibility of 
an internal compensating transformation. As a result, a 
violation of rotation symmetry for the metric correlation 
functions in the continuum becomes plausible, and we dis- 
cuss this issue in sect. IVIIII 

Furthermore, in our setting the existence of appropri- 
ate internal rotation or Lorentz transformations among the 
zweibein components e M ° and e^ 1 is not guaranteed a pri- 
ori. We investigate this question in the next section. There 
we find that no transformation of the field variables can 



account for a euclidean rotation among e M ° and e^ 1 . In 
consequence, we expect for the continuum limit the pres- 
ence of "preferred axes" and a violation of the continuous 
rotation symmetry. The issue of Lorentz symmetry is more 
subtle and will be discussed in the next section. 



VII. GENERALIZED LORENTZ 
TRANSFORMATIONS 



The formulation of the action (|47|) in terms of the 
zweibein suggests an investigation if a continuous Lorentz 
type symmetry could be present. Such a symmetry would 
have to act among the different components of the scalar 
field. Indeed, if we can find a transformation of the scalars 
such that (e^ ,^ 1 ) transform as a two-component vector 
(with index m), then the contraction with e mn in eq. (|41l) 
yields an invariant. This holds for "generalized Lorentz- 
transformations" corresponding to the groups 50(2) or 
50(1,1). (We employ here the name of "generalized 
Lorentz transformation" since they act on the zweibein in 
analogy to the usual Lorentz transformations. No fermions 
or spinor representations arc involved, however, in our set- 
ting.) 

Expressing the two complex scalars ifi in terms of four 
real components ip a , 



we can write the zweibein as 



with real antisymmetric 4x4 matrices 



0-10 
1 
10 



, a 



0-1 
-1 
10 
10 



(58) 



(59) 



0-1 0, 
Under infinitesimal transformations 

5if) = aA?P (60) 
the zweibein transforms as 

8e™ = 2Aaip(a m A + A T a m )d^. (61) 
A realization of the Lorentz group 50(1, 1) requires 

a A + A T <7° = a 1 , a 1 A + A T a 1 = o\ (62) 

such that 

5e^ = ae^ 1 , ^e^ 1 = ae M °. (63) 
The most general solution of eq. (f6"2"j) is 



.4 



(3 + (I- 7)r 3 , e + T]T 3 
C + X T 3 , - P + 7 T 3 



(64) 
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with arbitrary coefficients /3, 7, e, 77, £, x and Pauli matrices 
TV On the other hand, a realization of SO(2) would replace 
the first equation (|62[) by er'M+yl^V — —a 1 , while keeping 
the second one unchanged. No solution for A exists in 
this case. At this stage an internal Lorentz transformation 
remains possible, while internal euclidean rotations cannot 
be realized for our setting of scalar fields. 

For no choice of parameters /? . . . x the matrix A in eq. 
(|64|) is antisymmetric. The infinitesimal transformation 
(j6TJ|) does therefore not respect the condition 4> a ip a = 1 
for all allowed values tpa- It is n ot a genuine symmetry of 
our model. The basic reason is that 50(1, 1) is a noncom- 
pact group, while any constraint on the ip a that leads to 
a compact manifold (as the sphere S 3 in our case) admits 
no noncompact isometries. The issue may be understood 
in more detail by specializing to A = ^diag(l, — 1, 1, — 1) 
or 7 = 1/2, j3 = e = ry = C = x = 0. The transformation 
((60|) induces a change in the length of the vector ip a , 



(65) 



This does not vanish for arbitrary configurations which 
obey ip a ipa — 1- Nevertheless, it vanishes for the sub- 
class of configurations which obey ipf + ipl = ^>§ + V 1 ! or 
l^i 1 2 = \<fi2\ 2 - For example, this condition is obeyed for the 
stripe configuration (1231) . 

This observation has an interesting consequence. For 
any stripe configuration (|23j) . (|24j) we can infinitesimally 
increase all | ip\ | and decrease all | if2 \ , such that | <f\ \ 2 + \ if 2 \ 2 
remains unity for all lattice points. Keeping the phases (IM1) 
fixed this does not change the action. Such a transforma- 
tion amounts to an infinitesimal Lorentz rotation among 
the zweibein components e M m according to eq. (J63J) . 

We summarize that the action @ is invariant under 
Lorentz transformations of the group SO(l, 1) 



S(fi(x) = -a{x)<fx (x) , 6if 2 (x) 



1 



a(x)ip 2 (x)- (66) 



For the continuum action this symmetry is even a local 
symmetry, with transformation parameter a(x) depending 
arbitrarily on the position x. Indeed, the inhomogeneous 
part ~ d^a resulting from the derivatives cancels. The 
local character of the Lorentz symmetry is, however, not 
respected by the lattice regularization - the lattice action 
is invariant only for constant a. Furthermore, the non- 
linear constraint ((IJ is not compatible with the Lorentz 
symmetry. 



VIII. METRIC CORRELATION FUNCTIONS 
AND SYMMETRIES 

In this section we investigate the behaviour of fluctua- 
tions of the metric around the expectation value, 



h^(x) = g^u{x) - g^ v . 
The correlation function, 

G^ v ,pa{x,y) = (h f _ ll/ (x)h pa (y)), 



(67) 
(68) 



characterizes the response of the metric to a source (energy 
momentum tensor) in linear order. In our model we find 
in the metric sector a non-trivial scaling behavior without 
a finite correlation length. In Fig. [3] we exhibit the two 




FIG. 9: Euclidean metric correlator Gq 00 for different values 
of /3, parallel to the Eq axis. 
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FIG. 10: Diagonal euclidean metric correlators, defined in eq. 
(|69"ll. at P = 10. 



point function Goo, 00 for the euclidean metric (|16]) . We ob- 
serve a power law decay with the exponent a = 2 at larger 
couplings, while close to f3 c the power is slightly bigger. 

To disentangle this behavior we must decompose this 
correlator with respect to the discrete symmetries of the 
system. We list in Table U the behavior of the various met- 
ric components with respect to the discrete lattice symme- 
tries, with + for even and — for odd. Here the eigenstates of 
the reflection symmetries are denoted by goi, h = /loo + ^n 
and d = hoo = /loo — ^<5oo = 5(^00 — ^11 )• Correlation 
functions involving two powers of a given eigenstate are 
even with respect to all discrete symmetries. As result, the 
correlations 
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■|-rotations 


P,T 


D ± 


h = hoo + /in 
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+ 


d = §(/ioo - hn) 




+ 




hoi 
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(M) 


■|-rotations 
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+ 
+ 


+ 

+ 
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TABLE I: Symmetry properties of the euclidean metric with 
respect to the discrete symmetries of the lattice. 



TABLE II: Symmetry properties of the Minkowski metric with 
respect to the discrete symmetries of the lattice. 



Gh{x - y) 
G d {x - y) 



(h(x)h(y)) 

Goo, oo + Goo, 11 + Gn,oo 
(d(x)d(y)) 



Gi 



1 



= -AG 



00,00 



- G, 



00.11 



-G l 



1,00 



G ff(x-y) = (hoi(x)h 01 (y)) = G i, i, 



G1141), 
(69) 



are invariant if the difference y — x is rotated by an angle 
7r/2. Wc plot these correlation functions a function of dis- 
tance in Fig. [TO] We observe that the decay of Gd for large 
(x — y) is reasonably fitted by a power law 



Gd(x -y) = A d \x - y\~ 



(70) 



with a — 2, while Gh and G Q // show a short distance be- 
haviour with a steeper decay. The non-diagonal correlation 
functions of the euclidean metric described below show also 
a short distance behaviour The correlator Goo,oo(x,y) is 
thus a linear combination of the correlators above, some of 
which have short distance behaviour, and some has a long 
distance behaviour. For j3 close to /3 C the short distance 
contributions change the short distance behaviour slightly, 
while for larger j3 the correlator is dominated by Gd, which 
has a decay exponent independent of (3. 

In the continuum limit for r S> A the correlator Gd 
is the only non-vanishing correlation function on the Xq 
and Eq axes for the euclidean metric g)a> ■ The tensor 
properties of the metric imply that for a 7r/4-rotation Gd 
should be replaced by G // if rotation symmetry is realized. 
This is clearly not the case, therefore the euclidean rotation 
symmetry is broken in the continuum limit. The behavior 
of the correlation functions is related to an expansion of the 
effective action in second order in the fields . Our findings 
clearly indicate that the effective action is not of the simple 
form discussed in the Appendix [Al 

Using the discrete symmetries above one can find that 
certain correlation functions have to vanish, as described 
below. For example, the correlation function 

(d(x)h 01 (y)) = ±(G 000 i-G 1101 ). (71) 

is even under 7r/2-rotations. Since this correlation is odd 
under P, T and D± it has to vanish on the a; - and x 1 -axes, 



as well as on the diagonal axes 
hand, the correlations 



£0 or ~ Ei . On the other 



G h d{x -y) = (h(x)d(y)} = G 00 ,oo - Gi 



1,11 



(72) 



and 



Ghoff{x - y) = (h(x)g i(y)) = G 00 ,oi + Gn.oi (73) 

change sign if y — x is rotated by tt/2. The correlation G^d 
vanishes if y — x is parallel to one of the diagonal axes ~ Eq 
or ~ Ei , while Ghof f is zero if y — x is parallel to the a; - 
or a^-axis. These features are confirmed by our numerical 
results. 

Equivalently, we can investigate the action of the discrete 
symmetries directly for the correlations (|68|) . For example, 
the correlation functions of the type 



G^oiiy) = (hnv(0)hoi(y)} 



(74) 



have to obey 



GoooiCTy) = G ooi(Py) = -G ooi(y), 
Gnoi(Ty) = Guoi(Py) = -Guoi(y), 
Goioi(Ty) = Goioi(Py) = Goioi(y), 



(75) 



for 



y = (y°, y 1 ) , Ty = (-y°, y 1 ) , Py = -fa). (76) 

In particular, the correlations G0001 and Gnoi have to van- 
ish on the x°- and a^-axes, i.e. for y = (y°, 0) or y — (0, y 1 ). 
For y on one of the diagonal axes, e.g. for y = (m, m) or 
y = (m, —to) this implies for the correlations 



(h 00 (0)d(y)) = -(h n (0)d(y)), 
(h O i(0)d(y)) = 0. 



(77) 



The symmetry arguments are the same for all euclidean 
metrics - they hold both for h^J and . 

The correlation functions for the Minkowski metric be- 
have similarly to the euclidean correlations. The decompo- 
sition is now done whith respect to a Mikowski signature, 
as displayed in Table [TTJ In Fig. [11] we show the diagonal 
correlations 



G H {x - y) 
G D (x - y) 



(H M (x)H M (y)} 



G, 



(M) 



r 



(M) 



_ sty" 1 
00,00 ^00,11 ^11,00 

(D M (x)D M (y)) 
t(G 



r (M) 



AM) , r (M) 
T 00.00 ~r ^00,11 



^11,00 



^ll.llJ: 



Go f f(x-y) = (C ) (^)^i JJ (y)> = G r; , 1 



,(M) 



(78) 
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Similarly to the euclidean case, Gd shows a long dis- 
tance behaviour of a power law decay with a ~ 2. Along 
with GD,off this is the dominant correlation function in the 
continuum limit. The correlation function G // for the off 
diagonal metric components also shows a power law decay 
and remains present in the continuum limit, even though 
its value is smaller the Gd and Go,off by an order of mag- 
nitude. On the other hand, Gh decays quickly and plays 
no role in the continuum limit. In Fig. [12] we show the 
correlator Gq^q which behaves similarly to Gqq 00 shown 
in Fig. M 
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TABLE III: Symmetry properties of zweibein fluctuations. 
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FIG. 11: Diagonal Minkowski metric correlators, defined in eq. 
(75). for $ = 10. 




of j3, parallel to the Eq axis. 



IX. ZWEIBEIN CORRELATIONS 

The intriguing features of the correlations for the eu- 
clidean and Minkowski metric can be understood better in 
terms of the correlation functions for the zweibein fluctua- 
tions 



As we have discussed above the expectation value e^™ pre- 
serves the discrete symmetries P, T, D± and 7r/2-rotations. 
We classify the zweibein fluctuations according to these 
symmetries 

s = litf + fi 1 ) , a=^(/o°-/i 1 ), 

u = \{k l + h a ), v = ~(f x -f 1 °). (80) 

The transformation properties of the eigenstates s, a, u, v 
under P, T, D± and 7r/2-rotations are collected in tablc lllll 



This implies that the correlations (su) , (sv) , (au) and 
(av) must vanish on the a; - and x 1 - axes, while 
(so) , (sv) , (ail) , (uv) are zero on the diagonal axes ~ Eq 
or ~ Ei. A rotation of the axes by 7r/2 leads to a minus 
sign for the correlations (sa), (su), (av) and (uv), while the 
other combinations are invariant. 

The zweibein correlators also show a powerlaw decay 
similarly to the metric correlators, as shown in Fig. 1131 
We also plot the non-vanishing non-diagonal correlations 
on the Eo and xq axes in Fig. Q31 Again, only some modes 
exhibit the slow, power law decay, while other modes show 
faster decay, similar to the correlations of the metric ten- 
sor. The correlations relevant for the continuum limit are 
(aa), (vv) and (av). 
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FIG. 13: Diagonal correlations of zweibein fluctuations on the 
Eq axis measured at ft = 7. 
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FIG. 14: Non-vanishing non-diagonal correlations of zweibein 
fluctuations on the Eq and xq axes measured at = 7. 



FIG. 15: Correlations of the euclidean metric compared with 
estimates from zweibein correlations according to eq. (|86jl mea- 
sured at P = 10. Since Gq^' 01 has an alternating sign, we show 
the absolute value. 



The fluctuations of the metrics (1501) are related to the 
zweibein fluctuations by 

h(S) _ p rn~ n~ _ i~ m~ n\~ 



= {U m f~ 



m £ n _ £ ra n , m n 



(81) 



Taking the expectation value of eq. (1811) and employing 

(hffj) = 0, (ffj, m ) = 0, relates the difference between the 
metrics (fl~3|) and (l54| to the local zweibein fluctuations. 
This holds both for euclidean and Minkowski signature, 
such that for S = M, E one finds 



a (S) _ JSe) If mf n\fi 
■)p.v -lp.v \J/i Ju I 'In 



(82) 



Furthermore, if we neglect higher correlations ~ 

(/ 3 ),(/ 4 ) c one obtains an approximate relation between 

(s) 

the metric correlations for g^J and the zweibein correla- 
tions. (Here (J 4 ),- ((/ 2 (x) - {f\x))) (f(y) ~ (/ 2 (y)))>, 
with indices and fj omitted.) This leading order relation 
reads 



pp ^vm^ar 

4- F nr e p 

' vp ^pn^-ar 



Fva e p.n€ps, (83) 



where we define 



and 



rpmn 
pu 



-'fin 



(84) 



(85) 



In Fig. [TS] we show that the metric correlations are 
well aproximated using the estimate (|86|) in terms of the 
zweibein correlators. This is true even for the case of a 
fastly decaying mode. This evidence means that the higher 
order zweibein correlations quite small. It also implies that 
the power-law behaviour of the metric correlations is gov- 
erned by the power law behaviour of the zweibein corre- 
lators. In Fig. [IB] we show the (aa) correlator for differ- 
ent values of /3. Within our precision the decay exponent 
a ~ 2 is found to be independent of /3 in the stripe phase 
for/3>/3 c . 
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FIG. 16: Correlations (aa) of zweibein fluctuations on the Eq 
axis for different values of j3. 



For the euclidean metric correlation this yields in the 
striped phase given by eq. (|2"3")l 

r<( B ) — /ipoo ME) _ AP ii 

Wo, oo — ^00 ' "11,11 — ^ r ll i 
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4F, 
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4f; 



00 
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""01,01 ~~ ^11 



■2F- 



10 



F 



oo • 



(86) 



Finally, the scalar correlations are qualitatively different 
from the metric and zweibein correlators. Some typical 
scalar correlators are shown in Fig. [TTj As one observes, 
they can be approximated by an exponential decay rather 
than a power law. The decay constant is roughly indepen- 
dent of the coupling above /3 C , and the overall amplitude 
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of the correlations decays with increasing j3. For the scalar 
correlator we must take into account that the ground state 
in the broken phase is only invariant with respect to trans- 
lations by 4-\/2A. Therefore we plot the scalar correlator 
for distances which are a multiple of 4^/2A. 

We have shown in Fig. [7] that the metric expectation 
value is well approximated by disconnected contributions, 
i.e. the zweibein expectation values. In the first part of 
this section we have shown that also the metric correlators 
arc well described by the zweibein correlators, neglecting 
higher order zweibein correlators. In this sense the metric 
correlations are largely determined by the zweibein. The 
relation of the zweibein and the scalars seem to be differ- 
ent. The zweibein expectation value is well described by 
the disconnected part, i.e. the scalar expectation values, as 
shown in Fig. [5] One may ask if this can be extended to the 
correlation functions by neglecting higher order connected 
correlation functions of the scalars. We address this issue 
in appendix [C] However, the different qualitative behavior 
of zweibein correlators and scalar correlators make it un- 
likely that the zweibein correlator can be approximated by 
a linear combination of the scalar correlators. 
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FIG. 17: Typical scalar correlators along the Eo axis measured 
at /3 = 4.9. Note that the x axis has a linear scale. 



X. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the non-linear a- 
model given by eqs. ([1}, This model can be inter- 

preted as a classical statistical model. In this case j3 plays 
the role of an inverse temperature. This new class of non- 
linear cr-models shows several interesting features. For all 
temperatures below the critical temperature, /3 > /3 C , we 
find a universal critical behavior. Our model is therefore 
an example for self-tuned criticality, where critical behav- 
ior occurs without a tuning of parameters. This contrasts 
with the usual situation where critical behavior occurs only 
for a special choice of parameters, e.g. f3 = f3 c . The pres- 
ence of long range correlations for all j3 > j3 c shows cer- 
tain analogies with the Kosterlitz-Thouless phase transi- 



tion . Our model belongs, however, to a new univer- 
sality class that differs in important aspects from the uni- 
versality class characterizing the Kosterlitz-Thouless phase 
transition. It is not restricted to an abelian symmetry 
5*0(2) and shows a very different long distance behavior. 
It remains to be seen if some particular condensed matter 
system realizes this universality class. 

A characteristic feature of our model is the stripe phase, 
where the symmetries of discrete reflections and ir/2- 
rotations on the lattice are preserved only if they are com- 
bined with gauge transformations acting on the flavor in- 
dices of the scalar field. The stripes single out preferred 
axes for quantities that are not invariant under gauge trans- 
formations. The presence of a stripe phase is not mainly a 
property of the lattice formulation. We show in appendix 
[Dl that a stripe phase is also present in a continuum model 
with effective action similar to eq. (2). Such an effective 
action is a candidate for the description of the universality 
class, while we have not yet made any detailed comparison 
of its properties with our numerical results. Nevertheless, 
we have established that the effective action discussed in 
appendix [D] can describe the first order phase transition 
from the disordered to the stripe phase. 

The most prominent property of our non-linear cr-model 
is lattice diffeomorphism invariance. This implies general 
coordinate invariance or diffeomorphism symmetry in the 
continuum limit. Our main interest concerns the possi- 
ble relations to quantum gravity and we therefore focus on 
collective order parameters that play the role of a metric 
and a zweibein. This realizes geometrical features for our 
statistical model, close to the conceptual framework of ref. 
fl3j . We have computed expectation values as well as cor- 
relation functions for the collective fields using numerical 
simulations. 

Our main findings are mentioned in the introduction and 
we may concentrate here on two aspects: (1) non- vanishing 
expectation values of the metric for the vacuum state, de- 
scribing flat space with either Minkowski or euclidean sig- 
nature, (2) long range correlations for the metric fluctua- 
tions that decay with r~ 2 as a function of euclidean dis- 
tance r. 

As in usual gravity, the geometry of space or spacetime 
is deformed by the presence of matter. This can be seen by 
introducing an energy momentum tensor as a source term 
for the metric fluctuations. The response of geometry to 
a local energy momentum tensor is fixed in the linear ap- 
proximation by the metric correlation function. One finds 
for the deviation from flat space 

VW = / (87) 

Jy 

with t pa an appropriate energy momentum density and 
Gf_i»,po- the metric correlator defined in eq. ([SS]). In partic- 
ular, for t 00 (y) — KS(y) and vanishing other components, 
one obtains 

= KG„ u , oo (x,0). (88) 

For this source the metric perturbation hoo, that may 
be compared to the Newtonian potential in four dimen- 
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sional gravity, decays ~ r 2 , with euclidean distance r 2 = 
visible in Fig. |H1 
A static localized massive object in a Minkowski setting 
would correspond to a time independent energy momentum 
tensor 

t m (y°,y 1 ) = M6(y 1 ). (89) 

This results in a static metric that decays inversely pro- 
portional to the spatial distance \xi\, 

M , , 

hoo{x ,xi) ~ i — r, (90) 
Fi| 

similar to four dimensional gravity. 

Our model can be considered as a model for two- 
dimensional quantum gravity in the sense that the effective 
action for the metric is invariant under general coordinate 
transformations and that the metric correlations are long 
range. As a perhaps surprising effect the vacuum state 
corresponds to flat space even in the presence of quantum 
fluctuations. There are, however, also important differ- 
ences as compared to Einstein gravity in four dimensions. 
They are mainly related to the non-trivial stripe order for 
(3 > (3 C . The presence of preferred axes may lead to features 
that are not encountered if all order parameters preserve 
Lorentz symmetry as for standard four-dimensional grav- 
ity. The stripe order is reflected in the behavior of the 
correlation functions for the zweibein. In appendix E we 
discuss a simple ansatz for the zweibein effective action in 
the presence of stripes. It seems to differ substantially from 
the simple Lorentz-invariant setting of appendix B. While 
the effective action for scalars and zweibein is assumed to 
be Lorentz invariant, the Lorentz symmetry may be spon- 
taneously broken the vacuum by the stripe configuration. 

Several important issues remain to be solved before a 
more realistic model for quantum gravity can be con- 
structed from a suitable scalar field theory on a lattice. The 
approach using a collective vielbein seems quite promising. 
A diffeomorphism invariant lattice action can then easily 
be formulated by employing the determinant of the col- 
lective vielbein. In this case one would like to implement 
the Lorentz transformations acting on the internal or fla- 
vor index of the vielbein as an exact symmetry. This is 
possible along the lines discussed in appendix [D] One may 
employ a Lorentz invariant potential V{p\p2) already for 
the microscopic lattice action. The non-linear constraint 
([1]) can then be replaced by a Lorentz invariant constraint, 
for example by a bound on the Lorentz invariant product 
P1P2 < C. The second issue concerns the preservation of 
a global Lorentz symmetry for a flat space ground state. 
One possible solution is a non-zero expectation value for 
the vielbein e™ ~ <5™> while the order parameter for the 
stripe configuration vanishes. Finally, an important step is 
the transition from two to four dimensions. 

Several of the mentioned problems are absent or solved in 
lattice spinor gravity. However, reliable computations are 
difficult for spinor gravity. For this reason it seems worth- 
while to pursue in parallel the scalar approach to lattice 
gravity which permits relatively cheap numerical simula- 
tions. 



Appendix A: Propagating metric in two-dimensional 
gravity 

In this appendix we demonstrate that two-dimensional 
gravity can have propagating metric degrees of freedom. 
The issue which degrees of freedom propagate depends on 
the form of the quantum effective action. We present a 
simple example for such an action where metric degrees of 
freedom are indeed propagating. 

In two dimensions the curvature tensor has only one 
independent component that is related to the curvature 
scalar R by 

#0101 = g det(g^)R. (Al) 

The integral J y/gR is a topological invariant and does 
not contribute to the field equations of the metric. If only 
this term and a two-dimensional cosmological constant are 
present in the effective action there will be no kinetic term 
for the metric, such that the metric is not a propagating 
field. 

Nevertheless, a diffeomorphism invariant effective action 
can be constructed as 

T = / JgRf(-D 2 )R, (A2) 

J X 

where D 2 = D^D^ and is the covariant derivative. This 
can be generalized by adding terms with higher powers of 
R. We will concentrate on non-local invariants [14j of the 
type 

r = cJ K , I K = [ ^/gR{-D 2 )- K R, (A3) 

J X 

with n > 0. Our setting remains more general, however. 

In linear order of an expansion around flat space, /i M „ = 
9fif ~ ViJ.v-,h — h^yXf" 11 one finds (for arbitrary dimension 
d) 

R = d ll d v h^ -d 2 h = -^i-(9 2 C, (A4) 

with 

C = h - ^T^r^ - V = V + \*nw (A5) 

It is easy to check that £ is invariant under the inho- 
mogeneous part of the gauge transformations, S^h^ = 
dfj,£v+dv£ij,- For d = 2 the invariant (|A3|) reads in quadratic 
order in h^ v 

I K = \Jj(-d 2 f-«C (A6) 

For an effective action V = c K I K , K > 1, and for euclidean 
signature rj^ = 5^, this implies correlation functions de- 
coupling for large (x — y) as 

{C{x)ay))-\x~y\ 2 - 2 \ 
(d 2 ax)d 2 C(y))~\x-y\- {2+2K) . (A7) 
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For Minkowski signature the field equation for £, which is 
obtained by taking a funtional derivative of the effective 
action (|A6j) . describes a relativistic wave equation for a 
propagating degree of freedom. 

A non-local effective action of the type (|A3[) typically 
indicates the presence of a massless degree of freedom. For 
a suitable choice of degrees of freedom, T can often be 
written in an equivalent (quasi-)local form. For example, 
an effective action 

r = j^^-\xZ{-D 2 )(-D 2 )x + fxR^ (A8) 

describes a scalar field x with non-trivial wave function 
renormalization Z of the kinetic term and local coupling 
to the curvature scalar R. This implies for the scalar field 
X the field equation 

X = fZ-\-D 2 )(-D 2 )- 1 R. (A9) 

Insertion of eq. (|A9[) into eq. (|A8I) yields the gravitational 
effective action 

r = : y J x V9 Rz ^(- D2 )(- D2 y lR - ( Ai °) 

For Z = —D 2 we recover the action (|A3|) with k = 2. 
For an expansion around flat space the scalar x is directly 
related to £ in eq. (|A5|) 

x = f -Z-\-d 2 )C (All) 

The scalar x or C is the only propagating field in this type of 
two-dimensional gravity. Correlation functions for metric 
components with an overlap with this scalar field, e.g. for 
hoo or /in, should show a powerlike decay given by eq. 
■ fATl; . 

The effective action (|A3|) is invariant under general co- 
ordinate transformations. This is reflected in the linear 
expansion by the fact that only the particular combination 
C of metric components contributes to the action. Without 
diffeomorphism symmetry nothing particular distinguishes 
£ from the other components of h^v- The expansion of 
the effective action (IA3|) around a flat space exhibits the 
symmetry of global rotations or global Lorentz transfor- 
mations, depending on the signature of the metric. In this 
important aspect it differs from the effective action for the 
collective metric in our non-linear a- model 

Appendix B: Possible effective action for propagating 
zweibein in two dimensions 

In this appendix we briefly discuss an example for a pos- 
sible form of an effective action for the zweibein. This 
serves as an illustration of some of the effects that may 
be expected for a more realistic effective action. We in- 
sist on diffeomorphism symmetry of the effective action, 
but we do not require local Lorentz symmetry acting on 
the index m of the zweibein e™. Similarly, we also do not 



impose euclidean rotation symmetry. As a consequence, 
the covariant derivative contains no spin connection, and 
Z? p e™ can differ from zero. (A discussion of this type of 
generalized geometry can be found in ref. [HI). 

Our example for an effective action involves the deter- 
minant of a "renormalized zweibein" 

„Rm 71111/ n IT3-[\ 

e n - -Vt e ^ • ( m J 

Here the "wave function renormalization" ZJ^ is a func- 
tion of covariant derivatives D M such that e Rm transforms 
again as a covariant vector. An example is (D 2 = D^D^) 

Z™ = Zl (-D 2 )6;6™ + z 2 (-D 2 )D,D»S™, (B2) 

with z\ and z 2 scalar functions. (Covariant derivatives in- 
volve the Levi-Cevita connection in the usual way, but no 
spin connection.) For the leading term in the effective ac- 
tion our ansatz reads 

r — —e^f [ d 2 re Rm e Rn 

i 4A mn I n v 

with 

A% n = Z;% l ZZe p «e rs . (B4) 
For the example (IB2[) one has 

A%n = \A.^" + z lZ2 {D p D^P" -DpD'e^) 

+z 2 .D p D^D r7 D»e p °}e mn . (B5) 

The effective action (|B3j) is diffeomorphism invariant 
provided that covariant derivatives of tensors are again ten- 
sors, such that A^ n transforms as a scalar function multi- 
plied by & v . The effective action would be invariant under 
a global generalized Lorentz transformation if the covari- 
ant and contravariant derivatives are singlets with respect 
to this transformation. This is possible only for one of the 
groups 5*0(2) or 50(1, 1), but not for both simultaneously. 
For example, a definition 

D" = e m »e n v f, mn D v = g^D v (B6) 

requires the specification of f] mn , with inverse fj mn obeying 
V mn Vnp = S™. The contravariant derivative also involves 
the inverse zweibein e m M which is defined by the relations 

p n — t) n p m p u = S v (B7) 

In order to gain some intuition for the implications of the 
effective action (|B3[) we first consider the approximation 
where z\ is a function of — d 2 = —d ll d u f) ,lv , while 22 = 0. 
The quantum field equation reads 

^ = -§^-^ = 0. (B8) 
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Flat space with e„" ~ 5™ solves this field equation provided 
that the Fourier transform z\ (q) vanishes for vanishing mo- 
menta q = (qo,qi). The second functional derivative reads 
in momentum space 



S 2 T 



Se ll m (-q)Se^(q') 



M 

2A 2 ' 



r'"' 



S(q,q') 



(B9) 



For z\ a function of q 2 = q^q^ = —q 2 +q 2 the corresponding 
inverse propagator is a Lorentz covariant expression. 

In the space of the four-component vector E = 
(eo°, eo , ei , ei 1 ) the inverse propagator takes the form 
(we omit the (5-function in Fourier space which reflects 
translation symmetry) 



G~ 



M 

2A 2 



B. 



(BIO) 



with B an invertible 4x4 matrix 



B = 



fo 








1\ 








-1 








-1 



















£? = 1. 



(Bll) 



The propagator in momentum space is therefore given by 



G 



2A 2 



The non- vanishing correlation functions 
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(B12) 



sfV)- (B13) 



This clearly differs from the observed structure of zweibein 
correlations for our non-linear cr- model, as discussed in 
sect. HXl It becomes clear however, that very different struc- 
tures can also be accounted for by the ansatz (IB3|) . In ap- 
pendix [E] we discuss a possible form of an effective action 
for scalars and zweibein that may be somewhat closer to 
our model. 



namely 



^{ x ))^ 5 {y))^^A Hai {x,y). (C4) 



The relation (|C4[) between FJ^J 1 and H ai involves the 
expectation values of the scalar fields. We will evaluate 
them for the stripe configuration (|21l) . (f2"3"|) . Since we are 
interested in large separations (x — y) we employ the con- 
tinuum limit. In the complex formulation we take for the 
cell averages continuous fields 

(piOr)) = ^=(l- l )exp|-^(.T°+x 1 ) 

{Mx)) = 272 (1 ~ l) exp {~jk ( ~ x ° + xX) } - (C5) 

These fields are indeed invariant under the combined trans- 
lations t and ti, cf. eqs. (f25|) . (f26| . Taking partial deriva- 
tives of eq. (|C5j) reproduces the relations (|28|) up to a factor 
7r/4. This conversion factor for derivatives between the dis- 
crete and continuum formulation (e.g. discrete derivatives 
involving finite distances versus continuous derivatives, and 
cell averages versus continuous fields) has to be applied to 
eqs. (jCip . (IC2[) . (|C4[) if we use the standard partial deriva- 
tives, i.e. dfj, -t (4/7r)5 M . (The factor 4A 2 in eq. (|C4|) gets 
replaced by (64/7r 2 )A 2 and the factor 2A in eqs. (jCljl . 
(|C2|) becomes 8A/7T. With this replacement the evalua- 
tion of eq. (|Cip for the stripe configuration (|C5|) yields 
indeed e p m = S^ m .) We also observe the normalization 

(<pt)(pi) = (</**){<&) =1/4- 

We finally take into account that the normalization of 
the expectation value differs from eq. (JlJ by multiplying 
(<p a ) by a factor Z v , such that |(ipi)| 2 + \ (tp2)\ 2 = Z 2 < 1. 
In terms of the real fields (-0a) the stripe configuration 
becomes 



Appendix C: Zweibein correlations from scalar 
correlations 

For large (3 one expects that the zweibein is well approx- 
imated by the scalar expectation values (ip) , 

e^ m = 2A(^ ct )< /3 ^(^). (CI) 

Similarly, one may try to approximate the zweibein fluc- 
tuations by the scalar fluctuations, Sip a — ip a — (ipa), by 
linearizing eq. (|58p. 

U' n = e™ - er = 2A{<5V> a <^ M <^) + 

(C2) 

In this approximation the zweibein correlations are approx- 
imated by the scalar correlations 

H a p(x,y) = (5ip a (x)Sipf3(y)), (C3) 




Inserting eq. (|C6[) into eq. (|C4|) yields explicit expressions 
for the zweibein correlations as linear combinations of the 
scalar correlations. So far, we have not attempted to check 
this type of relations numerically. The different qualitative 
behavior of scalar and zweibein correlations sheds doubts 
on the validity of such an approximation. 
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Appendix D: Effective action for scalars and field 
equations 

The quantum effective action T[ip] for the scalar fields 
is defined in the usual way by introducing sources for the 
scalar fields in the functional integral, and performing a 
Legendre transform of the generating functional for the 
connected Greens functions. It includes all effects of fluc- 
tuations and generates the one-particle-irreducible Greens 
functions. Thus the functions for an arbitrary number of 
fields follow from T by simple functional differentiation. In 
this sense the knowledge of T amounts to a solution of the 
model. We do not attempt here a computation of the ef- 
fective action. We rather investigate a simple ansatz which 
respects the symmetries of our model, namely 



1 



d^VW^a^d^^sd^s, (Dl) 



where the cr^g matrices are defined in eq. (|59")l . (We use 
ip a instead of (tp a (x)} in the following.) Eq. (IDlj) equals 
the classical action ([2]) for V(ip) = —(3. However, we admit 
here a general "scalar potential" V(iji). If V depends only 
on pi = ip\ipi = ipf + ipl and p 2 = ^2^2 = Vi + 
with V(pi,p2) = V(p2,pi), the effective action shares all 
symmetries of the classical action. We use the continuum 
version of the effective action (|D1I) in order to demonstrate 
that the phase transition to the stripe phase also occurs in 
a continuum theory. 

The vacuum state (or thermal equilibrium state in case 
of a classical statistical interpretation) is a solution of the 
quantum field equations. In our case, the quantum field 
equations for the scalar fields are obtained from the func- 
tional derivative of the effective action (|Dlj) . In the absence 
of sources they read 

+ d,[V^ 5 d^ 5 ]< p M = 0- (D2) 

Here we use 

dV 



di\> a 



(D3) 



with V' a = dV/dpi for a = 1, 3, V£ = 0V/dp 2 for a = 2, 4. 
The field equations always admit the solution ip a — which 
corresponds to the disordered phase. We are interested 
here in the stripe solutions which are obtained for 

2Va^d^^ v e mn e v n = -2V^ a , (D4) 

with V, V' a and 

e™ = V 7 <A^ , e = det(e AI m ), (D5) 

independent of x. Multiplication of eq. (ID4[) by ip a and 
summing over a yields as a condition for the existence of 
this type of solution 



dV dV 
Pi^ + P2j— = -2V. 
dpi dp 2 



(D6) 



We make the ansatz 



(D7) 



where only the elements (1, 3), (3, 1), (2, 4) and (4, 2) of the 
constant matrices Uo and Ui differ from zero. Eq. (|D4I) is 
obeyed for tp a ^ if 



e u „e mn e n u (2VK> m + V^ m )=0. 



Here u^' m is defined by the condition 

rn [rj \ „,a,mr 
cr ce/3{ U n)0'y = V d Q7 



and obeys 



The squared matrices {U^) 2 are diagonal 

{U^ap — C^,a5 a /3 , Cp4 = Cm, 3 > Cm, 2 = Cm, 4- 

Our ansatz (|D7|) requires then 

(<9 M ) 2 V>a = Cm, "V'a- 



(D8) 
(D9) 

(DIO) 

(Dll) 
(D12) 



We will require the matrices to be antisymmetric, guar- 
anteeing d^pi = d^p2 = 0. Then the coefficients Cm, a are 
negative (or zero). As a consequence, one obtains solutions 
of eq. (|D12[) which are periodic in both x° and x . We typ- 
ically will find solutions with Co. a = Ci,a- They obey the 
wave equation 



{dl - %)1> a = 0. 



(D13) 



It is remarkable how wave equations with two-dimensional 
Lorentz symmetry arise in a natural way from the field 
equations derived from the action (|D1I) . 

In analogy with eq. (IC6|) we consider possible solutions 
of the type 



^3 
V>2 
1p4 



with pi 
obeying 



r y cos (Pi (a; + x 1 )) - sin (Pi(x° + x 1 )) 

-ci sin (Pi(x° + x 1 )) + cos (Pi (a; + a; 1 )) 

c 2 [cos (P 2 (-x° + x 1 )) - sin (P 2 (-x° + x 1 )) 

-C2 [sin {P 2 {~x Q + x 1 )) + cos (P 2 (-x° + x 1 ))] , 

(D14) 

2c\, p2 = 2c 2 and antisymmetric matrices 

(U )l3 = Pi , (tfl)l3 = P l- 
(Uo)24 = —P2 j (C^l)24 = P2, 



Co, 



Ci,i = -Pi , Co,2 = Ci,2 = -Pi (D15) 



For a = 1,3 one finds for all p and m that u" ,m = Pi 

whereas for a = 2,4 one h as u" ' m = Pa if fi = m, and 
u a,m = _p 2 jf ^ ^ m Eq (jpXoj) yields 



PlPl + P2P2 , PlPl - Pj/92 

P1P1 - P2P2 , P1P1 + P2P2 



(D16) 
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For the particular stripe solution (|C6|) with Pi = P A 
7r/(4A) , pi = p 2 = one recovers 



W(p ) 



_ f_ rm 

~ 8A " 



OA 

e m — p m = Z 2 8 m 



For general Pi,P 2 ,pi, pi one has 



e = det(e M m 
and eq. (|D8|) is obeyed for 



4P 2 ^2PlP2, 



dpi dp 2 



(D17) 



(D18) 



(D19) 



If eq. (|D19[) has a solution for suitable values of p\ and pi 
we therefore find solutions with arbitrary Pi and Pi. The 
condition (|D19[) implies the condition (|D6[) . We conclude 
that for a potential which admits a solution of eq. (|D19[) 
stripe solutions (|D14[) exsits with arbitrary "momenta" Pi 
and Pi. 

We may try to interprete the effective action (|D1[) as 
an approximation to the continuum limit of the quantum 
effective action which corresponds to the microscopic lat- 
tice action @. In this case one expects a dependence of 
the shape of V on the parameter j3. (This extends to a 
parameter dependence of V for other models in the same 
universality class.) A phase transition from the disordered 
phase with ip a = to the stripe phase occurs at /3 C if for 
(3 > f3 c eq. (|D19[) has a solution and if for the correspond- 
ing stripe solution (ID14I) the action (|D1|) becomes negative. 
At the phase transition for (3 = f3 c the effective action in 
the stripe phase vanishes, such that the free energy T has 
the same value for the disordered and the stripe phase. A 
first order transition is realized if for /3 = /3 C the stripe 
solution still has a nonvanishing "order parameter" ip a . 

In the remainder of this appendix we discuss simple 
shapes of the potential V that realize the first order phase 
transition that we observe in our numerical results. A con- 
stant potential, e.g. V — — /3, does not admit solutions 
with ip a ^ 0. This has a simple explanation: the action 
is then a pure quartic polynomial of ip, such that for any 
value of ip for which T < the rescaled field (1 + e)^, e > 0, 
leads to an even smaller value of T, thus excluding an ex- 
tremum for T ^ 0. For the microscopic action this problem 
is cured by the non-linear constraint (JT|), which would be 
translated to the continuum language as p± + pi — 1/2. 

Solutions of the condition (|D19I) exist for a wide class 
of non-trivial potentials V(pi,pi) without invoking con- 
straints for pi and pi. As a first example we consider a 
potential V(p), p = pi + pi. The stripe solutions corre- 
spond then to an extremum of the combination 



W(p) 



P 2 V( P ) , f^(po) 
dp 



0. 



(D20) 



Indeed, eq. (ID20I) implies that the condition (|D19[) has a 
solution with P\ = pi = Po/2. (For j3 — > oo one would 
expect po — > 1/2.) The value of the effective action for 
stripe solutions with pi = pi = po/2 is given by Wq = 



Tn = 



J d 2 xeV(p ) = PiPi J d 2 xW . (D21) 



For P\ = Pi = 7r/(4A) the action per lattice point equals 
(tt 2 /8)Wq, such that for (3 —> oo one expects Wq —> 
— (8/n 2 )(3. We observe that Tq in eq. (|D21j) can be made 
arbitrarily negative for P1P2 — > 00. This "ultraviolet di- 
vergence" is cut off by the lattice regularization. We may 
consider (tt/AA) as the maximal momentum, say in the 
a; 1 -direction. 

Many different forms of the potential V(pi,pi) are con- 
ceivable. For example, if V only depends on the combina- 
tion pi pi the effective action (jDll) is invariant under local 
Lorentz transformations 

V4 ;3 (z) = e a ^^ h3 (x) , i>' 2A (x) = e-^/ 2 ^ 2A {x) 

(D22) 

for which the zweibein transforms as 



(e°)' 
(e 1 )' 



cosh a e 
cosh a e 



sinh a e 
sinh a e 



For the example (with positive constants a, b) 
V = -a + bpipi 



(D23) 



(D24) 



one has pidVdpi = pidV/dpi = bpipi and eq. (|D19I) is 
met for 



with 



Pi, 0^2,0 



V{Pw,P2o) 



2b' 



a 
2' 



(D25) 



(D26) 



Indeed, the combination eV takes for the configurations 
(jDT4)) the form 



eV = 4:PiP 2 [-apipi + b{pip 2 y 



(D27) 



For P1P2 > this has a minimum for eq. (|D25|) with e > 
0, (eV)o = — PiP 2 a 3 /(2b 2 ). We observe a degeneracy of 
the minimum under global Lorentz transformations (|D22[) . 

For P1P2 < the combination eV has a maximum with 
e < 0, eVb > 0. The action can become arbitrarily negative 
for P1P2 < and large pipi- We may prevent this to 
happen by imposing a Lorentz invariant constraint 



(D28) 



a 

P1P2 < 



An interesting issue concerns the possibility to use eq. (ID1I) 
with eq. (|D24[) for the microscopic action S, and to replace 
the non-linear constraint (p} by the condition (|D28[) . This 
would permit to realize Lorentz symmetry as an exact sym- 
metry of the model. 

For our model the observed first order transition at j3 c 
indicates an effective potential that is more complicated 
than the form (|D24p . Indeed, a first order transition is 
described if a possible local minimum of eV with pi, pi 
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different from zero occurs for j3 c at (eV^o = 0, while for 
(3 < f3 c one has {eV) > 0. For (3 < (3 C the absolute 
minimum of the effective action (|D1|) will then be found at 
ip = 0, corresponding to the disordered phase. 

As an example, consider 



c d 
V = -a + 6pip 2 - -(pi + p 2 ) + + p 2 ) 2 . (D29) 



Possible solutions of eq. (|D19[) with pi = p 2 = p/2 occur 
for 



with 



3c±^ + 32a(d + b) 

Po = WTb] (D30) 



V = V(p )=~-^cp . (D31) 



For d + b > and c > one finds a critical value 



4(d + 6)' 



(D32) 



such that for a > a c one has Vo < 0, and for a < a c the 
solution (|D30|) with smallest V occurs for Vo > 0. Thus for 
a < a c the disordered phase with tf> — is realized and we 
may associate the critical value /3 C with a(/3 c ) = a c . Indeed, 
for p\ = p2 = p/2 we can consider 



W(p) = p 2 V(p) = -ap 2 - C -p* + i(d + fe)p 4 . (D33) 



For a < one has a local minimum of at p = 0. 

A second local minimum exists for c > (4/3)y— 2a(d + 6). 
For c > 2yJ —a(d + b), corresponding to a > a c , this second 
minimum occurs for negative W and is therefore deeper 
than the minimum at p = 0. At the critical value a c the 
order parameter ~ p jumps from zero to 



Po,, 



2c 
d + b' 



(D34) 



Having found a satisfactory description of the first order 
phase transition from the disordered phase to the stripe 
phase one may ask if a suitable shape of V can also account 
for the scalar correlation functions in the continuum limit. 
In principle, the inverse scalar correlation functions can be 
obtained from the second variation of the effective action, 
evaluated for the appropriate solution of the field equation 
(|D2|) . We have not yet performed a computation of the 
correlation functions that correspond to the effective action 
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Appendix E: Effective action for zweibein and scalars 

The effective action for scalars and zweibein can be de- 
fined by introducing appropriate sources 

Wfat] = In J ^exp{-S + ]TC(y)e7(y) 

y 

z 

dW 



SC(y) 
dW 



(El) 



drj a (z) 

and performing a Legendre transform 

r^,e] = -W + Y,trn»(y)e™(y) + Y<VJS)MS)- (E2) 

y 2 

This yields the exact quantum field equations 



difj a (z 



-= Va {z). (E3) 



The scalar effective action T[ip] discussed in the preceding 
section is obtained for i m M (y) = 0. It can be inferred from 
r[i/;,e] by solving the field equation <9r/<9e M m (y) = with 

solution e^™ 1 {y)[i/j] being a functional of ip. Then T[iJj] — 

rty >e c°M]. 

Due to lattice diffeomorphism invariance of the action 
the continuum limit of the effective action is invariant un- 
der general coordinate transformations @. Besides diffeo- 
morphism symmetry the effective action also preserves the 
discrete reflection symmetries of the lattice action as well 
as the continuous flavor symmetry. 

Let us try an ansatz for the continuum limit of the ef- 
fective action which is consistent with the symmetries 



1 



p - mn 

2 



d 2 x{ip a a^ l3 d^ip [Vi (ip)^a^ s d v ip s 

+V 2 (ip)e 1/ n ]+V 3 (iP)e ll m e™}. (E4) 

The field equation for the zweibein reads 
1 



(E5) 



For Va,3 ^ 0, and in the absence of sources i m M = 0, the 
solution is 



Y2_ 

2V 3 



ip 7 a™ s d v ips- 



(E6) 



The proportionality between e™ and ip 1 a" 6 d u ^$ is realized 
in our model for f3 > /3 C , as can be seen in Fig. [SJ The 
proportionality factor /N^ e ^ corresponds to — V2/2V3, 
evaluated for the appropriate p\P2 (cf. appendix [D|l . and 
taking the proper normalization of partial derivatives into 
account. 
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Insertion of eq. (|E6|) into eq. (|E4[) yields eq. (|D1[) with 

The field equation for the scalars (|E3|) for ?? Q = yield, 
after insertion of the solution (|E6[) . precisely eq. (|D2[) with 
V given by eq. (|E7|) . We can therefore take over the dis- 
cussion of the preceding section. In particular, for (3 > f3 c 



one finds the wave solutions (|D14|) . with 

ip a <pd^p ~ 5™ , e M m ~ <5™. (E8) 

Eq. (IE6I) explains our finding e^ m ~ <5™ for all /3 > 
/3 C , and that e^" 1 as well as ip vanish simultaneously for 
(3 < (3 C . We realize a flat space geometry without tuning 
of parameters for all (3 > f3 c . 
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